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First preliminary measurement of the direct CP- violating parameter A g by the NA48/2 ex- 
periment at CERN SPS is presented. Using more than 1.6 billions of charged kaon decays into 
three charged pions, the charge asymmetry in the K — » 7r ± 7r + 7r _ Dalitz plot slope, A g , has 
been measured to A g = (0.5 ± 3.8) x 10~ 4 . This result is more than an order of magnitude 
more precise than results of previous experiments. 



Introduction 



Violation of the CP symmetry, and especially direct CP violation in the decay amplitude, due 
to its subtle nature, is an important window into physics beyond Standard Model (SM). 

It took more than three decades since the discovery of CP violation in the neutral kaons by 
Christenson, Cronin, Fitch and Turlay™, until direct CP violation was definitively established. 
After an unconfirmed indication by NA31 ^, KTeV and NA48 have demonstrated in the late 
90's with high significance that direct CP violation exists in the decays of neutral kaons into 
two pions"™. In the year 2001, B-factory experiments Babar and Belle have found CP violation 
in the system of neutral B mesons™ and last year also the direct CP violation in B-decays has 
been demonstrated™. 

In order to explore possible non-SM enhancements to heavy-quark loops which are at the 
core of direct CP-violating processes, all manifestations of direct CP violation need to be ex- 
perimentally studied and measured. In kaons, besides the already measured parameter e' in 
Kl — > 7T7T decays, the most promising complementary observables are decay rates of GIM sup- 



°Talk given at XXXXth Rencontres de Moriond, ElectroWeak Interactions and Unified Theories, La Thuile, 
March 2005. 

6 On behalf of the NA48/2 collaboration. 



pressed rare kaon decays^ which proceed through flavor-changing neutral currents and the 
asymmetry between K + and K~ decays into three pions. 

The — > 3ir matrix element is usually described in a polynomial expansion of two Dalitz 
variables u and v. 

\M(u,v)\ 2 oc 1 + gu + hu 2 + kv 2 + 0(u 3 ,v 3 ) (1) 

where \h\, \k\ <C \g\ and 

S3 - S , S 2 ~ Si 

u = 5 — and v = = — (2) 

with Si = (px — Piri) 2 and sq = J2 s «/3 (i = 1, 2, 3). Index i = 3 stands here for the odd pion c . 
Slope parameters g can differ between the K + and only due to direct CP violation d . SM 
predictions for the corresponding asymmetry 

A 9 = 4^ (3) 
9 9 + + g~ 

vary between few 10~ 6 to few 10~ 5 ^. An analogous asymmetry of integrated decay rates, Ay, is 
expected to be more than an order of magnitude smaller. Several experiments^ have searched 
for the asymmetry A g . The precision reached in both vr^vr+vr - and 7r ± 7r°7r° decay modes so far 
is at the level of few 1CP 3 . Existing theoretical calculations involving processes beyond SM^^ 
predict substantial enhancements of the asymmetry A g partially within reach of the NA48/2 
experiment. 

NA48/2 is an extension of the original experimental program of the experiment NA48 at 
the CERN SPS which has successfully accomplished the main goal to establish and measure 
the direct CP violation in the decays of neutral kaons into two pions ^. The primary aim of 
the NA48/2 extension j s t Q measure the parameter A g in tt^-k^it^ and vr^vr -^ modes to 
~ 2 - 4 x 10" 4 . 

NA48/2 took data in two periods. In total, during the 50 day run in the year 2003 and 
the 60 day run in the year 2004 about 4 billions of K ± — > vr^vr+vr and about 200 millions of 

— > 7r ± 7r°7r° decays have been collected and taped. The total recorded data volume amounts 
to about 200 TB. This paper describes the analysis and the preliminary result based on about 
1.6 billions of — > 7r ± 7r + 7r _ decays taken during the first 2003 run. 



2 Description of the Experiment 

NA48/2 deploys a novel system of two simultaneous charged beams with opposite charges. This 
allows to record decays of K + and K~ at the same time achieving a significant cancellation of 
systematic effects in the charge asymmetry measurement. 

The beams (Fig.^) of charged particles are derived from protons from the SPS which impinge 
with an intensity of about 7 x 10 11 protons per 5 s pulse on a beryllium target of 2 mm diameter 
and 40 cm length at a zero angle of incidence. The central momentum of 60 GeV and the mo- 
mentum bite of ±3 GeV are then selected symmetrically for positively and negatively charged 
particles in the first achromat unit which splits the two beams in the vertical plane and recom- 
bines them on the same axis. Then both beams pass through a series of focusing quadrupoles 
and are again split and cleaned in the second achromat. The second achromat together with 
the KABES detector^3 serves also as a beam spectrometer. The K + / K~ production ratio is 
about 1.8 (the analysis is independent of this ratio). 

After the second achromat both beams follow the same path. After passing the cleaning 
and the final collimators they traverse the entire ~114 m long decay volume superimposed 

c Other two, even, pions have equal charges. 

d Due to absence of mixing only direct CP violation is possible in decays of charged kaons. 




Figure 1: Lateral view of the NA48/2 beam and detector. The vertical scale is strongly enhanced. 



with a precision of about a millimeter. This superposition symmetrises the acceptances and 
contributes to the reduction of systematic biases. The decay region is comprised in a vacuum 
tank. The whole detector apparatus downstream the decay volume is almost identical to the 
NA48 experimental setup. 

The detector used for the reconstruction of K — ► 7r ± 7r + 7r~ decays is the magnetic spec- 
trometer. The spectrometer is housed in a tank filled with helium gas at atmospheric pressure 
separated from the vacuum tank by a KEVLAR window. A thin evacuated beam tube, travers- 
ing the centre of the detector, allows the undecayed beam particles to continue in vacuum. 
Two drift chambers are located before and two after the central dipole magnet which induces 
a horizontal transverse momentum kick of about 120 MeV/c to all charged particles. The drift 
chambers have an octagonal shape with an area of about 4.5 m 2 . Each is made of four sets of two 
staggered planes of sense wires oriented along four 45° directions. The momentum resolution of 
the magnetic spectrometer is cr(p)/p = 1.0% © 0.044%p with p in GeV/c units. 

The K — > 7T zil 7t + 'K~ decays are triggered with a two-level trigger system. At the first 
level, the rate is reduced to few hundreds of kHz by requiring at least two hits in a scintillator 
hodoscope placed behind the magnetic spectrometer. The second level trigger, which consists of 
hardware coordinate builders and a farm of asynchronous microprocessors, reconstructs tracks 
using data from the drift chambers. At least two tracks are required to converge within 5 cm 
in the decay volume. In the majority of the analysed data, level one triggers rejected by this 
condition are examined further and accepted even if only one track is reconstructed by the 
second level trigger. In this case it is required that this track is kinematically incompatible with 
a -k^it® decay, assuming the decaying kaon had momentum of 60 GeV/c and was moving along 
the beam axis. The resulting trigger rate is about 10 kHz. 

The description of other components of the NA48 detector apparatus can be found else- 
where^. 



3 Measurement Method 

The measurement is based on comparing the u-distributions of K + and K~ decays. In case 
of the tt^tt^tt^ final state, given the actual value ^1 Q f g = (-0.2154 ± 0.0035), the ratio 
Nk+{u)/Nk-(u) is proportional with sufficient precision to (1 + Agu). A g = Ag/2g is extracted 
from a linear fit to the ratio Nk+{u)/Nk~(u)- 



The presence of magnets both in the beam sector (achromats, focusing quadrupoles, etc.) 
and in magnetic spectrometer introduces an unavoidable charge asymmetry of the apparatus. In 
order to equalise local effects on K + and K~ beams the achromat and quadrupole polarities were 
reversed during data taking on an approximately weekly basis. The polarity of the spectrometer 
magnet has been reversed every day e . The whole approximately two- week cycle represents a 
super-sample which is treated in the analysis as an independent data unit. In the period of 
2003, four super-samples have been collected. 

Each super-sample contains four K + — > h + -k + it~ and four K~ — > ir~-K~iT + samples with dif- 
ferent combination of achromat and spectrometer magnet polarities. The ratio R(u) is obtained 
as a product of four Nk+(u) /Nk-{u) ratios: 

R(u) = RusRujRdsRdj « R{1 + 4A 5 n) (4) 

where U represents a configuration in which K + beam runs through the upper beam path in 
the achromats and D the lower. The index S represents spectrometer magnet polarity in which 
decay products having the same charge as the corresponding beam are deflected to the right with 
respect to the direction of the beam (towards the Saleve mountain) and J to the left (towards 
the Jura mountain). A linear fit to Eq. results in two parameters, normalisation R and Ag 
from which A g is extracted. 

The quadruple ratio technique in Eq. Incompletes the procedure of magnet polarity reversal. 
It allows a three- fold cancellation of systematic biases: 

• beam line local biases cancel between K + and K~ samples in which the beam follows the 
same path; 

• local detector biases cancel between K + and K~ samples deflected toward the same parts 
of the detector; 

• as a consequence of simultaneous beams, global time- variable biases cancel between K + 
and K~ samples. 

This method is independent on the relative size of the samples with different magnet config- 
urations. On the other hand, the statistical uncertainty depends on the statistically weakest 
of the eight samples involved in Eq. 21 Further reduction of systematic biases especially due 
to presence of stray permanent magnetic fields (earth's field, vacuum tank magnetisation) is 
obtained by maintaining azimuthal symmetry in the acceptance. 

Using the method described in this section, the result remains sensitive only to time vari- 
ation of asymmetries in experimental conditions which have a characteristic time smaller than 
corresponding field-alternation period. 

Due to superposition of the two beams the measurement doe s n ot need a Monte Carlo 
calculation of the acceptance. Nevertheless, detailed GEANT-based^*^ Monte Carlo simulation 
has been developed as a tool for systematic studies. The Monte Carlo simulation includes 
full detector geometry and material description, simulation of time-variable local drift chamber 
inefficiencies and time-variations of the beam geometry and drift chamber alignment. 

4 Data Analysis 

Several stages of compaction and filtering were necessary in order to reduce the data to a size 
suitable for the analysis. At least three reconstructed tracks in magnetic spectrometer, loose 



e During the data taking in the year 2004 the period of the spectrometer magnet reversal has been decreased 
to few hours. 



acceptance and quality cuts as well as at least one good reconstructed three-track vertex are 
required in the pre-selection phase. 

Tracks are reconstructed from hits in drift chambers using the measured magnetic field 
map rescaled according to the recorded current in the analysing magnet of the spectrometer. 
Chambers and individual wires are aligned using data collected in special runs in which muons 
were recorded with spectrometer magnet off. 

The three-track vertex is reconstructed using track segments from the first half of the spec- 
trometer correcting the extrapolations for the small magnetic fields due to magnetisation of the 
vacuum tank and from the earth's field. Using measured momenta and track directions obtained 
at the vertex the invariant mass of three pions is calculated. The stray-field correction, calcu- 
lated based on a three-dimensional field map measured in the entire vacuum tank, reduces the 
azimuthal variation of the reconstructed invariant mass of ~1 MeV/c 2 by an order of magnitude. 

The invariant mass resolution is about 1.7 MeV/c 2 . Since the decay — > 7r ± 7r + 7r~ is the 
dominant three-track decay, the sample is background free. The tails of the invariant mass dis- 
tribution are dominated by events in which one of the three pions decayed and the spectrometer 
reconstructs the track of the resulting muonA Using Monte Carlo simulations, the tails due 
to pion decays have been shown to be highly symmetric between K + and K~ samples. Only 
far tails are reject ed b y the cut \m n7T7T — vtlk\ < 9 MeV/c 2 , where mx is the PDG value of the 
charged kaon mass^l The systematic uncertainty due to pion decays, limited by the statistical 
precision of the generated Monte Carlo samples, is 5(Ag) = 0.4 x 1CU 4 . 

The most important feature determining the acceptance is the beam tube traversing the 
centre of all drift chambers. In order to securely exclude the central insensitive areas, all three 
tracks are required to cross the first drift chamber at least 11.5 cm from the beam centre 
and the last drift chamber at least 13.5 cm. The latter cut takes into account the additional 
beam deflection of about 2 cm with respect to the centre of the last drift chamber due to the 
spectrometer magnet. This important acceptance cut is related to the beam centre rather than 
to the centre of the detector. The reason is, that the beam optics can control the mean beam 
position only to about ±1 mm. The actual beam position is continuously monitored to a much 
better precision by calculating the momentum-weighted centre of gravity of three pions at first 
and at the last drift chamber planes, independently for K + and K~ . In addition to the time 
variation of the beam position, also the dependence of the beam position on the kaon momentum 
(~ ±1 mm in horizontal and ~ ±1 cm in vertical direction) is taken into account. In this way 
the K + and K~ acceptances cancel entirely and no Monte Carlo calculation is needed to correct 
for their difference. A conservative limit on residual systematic uncertainty, 5(Ag) = 0.5 x 10~ 4 , 
was determined by studying the sensitivity to various acceptance definitions. 

The measurement of the pion momenta is based on the knowledge of the magnetic field in 
the spectrometer magnet and on the tracking information from the drift chambers. The relative 
variations of the current in the magnet can be controlled down to about 5 x 10~ 4 . Smaller 
variations are continuously corrected by forcing the mean reconstructed tttttt mass to the PDG 
kaon massif with relative precision of about 10 -5 . This is done by scaling the measured track 
momenta symmetrically for positively and negatively charged tracks. As this effect is charge 
symmetric, by collecting K + and K~ simultaneously it cancels in the ratio R(u). 

A difference between the reconstructed tttttt invariant masses of K + and K~ is an unam- 
biguous measure of a residual horizontal drift chamber misalignment. An uncorrected horizontal 
shift of the chambers can lead to a charge-antisymmetric mis- measurement of the momenta. Ob- 
serving this invariant-mass difference as a function of time, revealed significant, up to about 0.2 
mm, movements of the drift chambers between individual alignment runs. Correcting the mea- 
sured momenta by p' = p(l + qf3p), where q is the sign of the track charge and (3 is proportional 



In order to avoid charge-dependent biases, the muon tracks are not vetoed by the NA48 muon-veto system. 



Table 1: Preliminary limits on systematic and trigger uncertainties on Ag in units of 10 







Snprtrmnptpr alip'Timprit, 


0.1 


Analysing magnet 


0.1 


Pion decay 


0.4 


Calculation of u and fitting 


0.5 


Pile-up 


0.3 


Total systematic uncertainty 


0.9 


Trigger efficiency - level f 


0.4 


Trigger efficiency - level 2 


0.8 


Total trigger uncertainty 


0.9 



to the difference of the invariant mass between K + and K~ 9 , reduces the effect to a practically 
negligible level' 1 . 

A potential source of systematic bias is the trigger. Inefficiencies of different trigger com- 
ponents are studied and measured using control samples from low bias triggers collected along 
with the main triggers. The rate-dependent parts of inefficiencies are assumed to be charge sym- 
metric (simultaneous beams). The inefficiency of the first level was measured to be small, about 
7 x 10~ 4 and stable in time. No correction is applied and an uncertainty of 5(Ag) = 0.4 x 10~ 4 , 
limited by the statistics of the control sample, has been attributed to this part. The rate inde- 
pendent inefficiency of second level trigger varies with time between 0.2 and 1.8%. The main 
source of this inefficiency are local drift chamber inefficiencies which are more important in the 
trigger than in the reconstruction due to reduced redundancy. All samples are corrected by the 
measured second level trigger efficiencies as a function of u. The uncertainty is fully dominated 
by the statistics of the control samples. 

Further sources of systematic effects were studied and evaluated. Bias due to resolution and 
the calculation of u was studied by changing the way u-variable is calculated from measured 
track momenta. Using odd or even pion tracks changes the resolution as a function of the 
position in the Dalitz plot (u, v). Additional studies were performed excluding various parts of 
the Dalitz plot. The preliminary result quoted in this paper is obtained from the fit restricted 
in the u-interval between —1 and +1. Effects due to pile-up of signals from different kaons due 
to high instantaneous intensity were studied by comparing samples with different amounts of 
recorded extra tracks close in time to the event. Charge-asymmetric material effects have been 
found negligible by studying amount and composition of the material in front of and in the 
chambers and taking into account pion spectra. Track charge misidentification was evaluated 
from events with three tracks of equal charge. Table ^ summarises all systematic uncertainties 
attributed to the result. 

5 Preliminary Result and Outlook 

The preliminary result presented in this paper is obtained from the full sample of reconstructed 
— > ir^ir^ir decays collected by the experiment NA48/2 during the data taking period of 
the year 2003. Three independent analyses, which agree within uncorrelated uncertainties, have 
been performed and averaged. The result is calculated separately for each of the four super- 
samples and then combined taking into account correlated systematic uncertainties (Table [2J . 



9 A shift of last drift chamber by 1 ^m corresponds to about 1.5 keV change in the invariant mass. 
h In the 2004 data-taking the alignment runs were performed on a more frequent, weekly, basis. 



Table 2: Measurement of Ag in units of 10 -4 in individual super-samples and the corresponding statistics. Only 
uncorrelated uncertainties, statistical uncertainty combined with the uncertainty from the level 2 trigger efficiency, 

are quoted. 
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Figure 2: Left: Measurement of Ag in four super-samples. Middle, Right: Detector resp. beam-line asymmetries 
which cancel in the Eq.^Jand their comparison to the Monte Carlo simulation. 



The measurement stability as a function of super-sample is shown in the left part of the 
Fig. [3 All four measurements are compatible with each other with x 2 /ndf = 3.2/3. As a 
systematic check, in the middle part of this figure, the quadruple ratio in Eq. 0] is rearranged 
such that instead of four K + /K~ ratios, four ratios of samples in which even pions are deflected 
to the right are divided by samples with even pions deflected to the left in the spectrometer 
magnet. In this case, the physical quantity Ag cancels and the result is expected to be equal 
to zero in the absence of any residual left-right detector asymmetries. Similarly, the right part 
of the Fig. [21 reflects the asymmetry of the two beam paths. These control asymmetries, which 
cancel at first order in Eq. EI show that the cancellation of systematic biases due to residual 
time variable imperfections in the apparatus is at the level of few 10~ 4 and therefore second 
order effects are negligible. Moreover, the comparison with Monte Carlo simulations shows that 
all these apparatus asymmetries are understood in terms of local inefficiencies and beam optics 
variations. 

The combined preliminary result from all four super-samples is 



Ag = (-0.2 ± 1.0 stat . ± 0.%^.^,) ± 0.9 sys t.) x 10" 4 



(5) 



Converted to the asymmetry using the PDG value of the Dalitz slope g< 



(0.5 ± 2A st at. ± stat. {trig.) ± 2.1 syst .) X 10 

(0.5 ±3.8) x 10 -4 



-4 



(6) 
(7) 



This result is compatible with no CP violation and with Standard Model predictions but has 
more than an order of magnitude better precision than similar previous measurements^. Further 
improvements are expected in future. The uncertainty due to the trigger can be significantly 
improved by deeper study of inefficiencies as well as by emulating the trigger in Monte Carlo 
simulations. Also some of the systematic uncertainties can be improved by more detailed studies 
of involved effects and by using data collected in the period of 2004. These data are being 
analysed and will significantly reduce the statistical uncertainty by more than doubling the 
total sample. 

NA48/2 has collected also about 200 millions of K — > 7r ± 7r°7r° decays. This decay mode 
is disfavored statistically due to its lower branching ratio and lower acceptance. On the other 
hand, this disadvantage is practically completely compensated by more favorable population of 
the Dalitz plot leading to an expected statistical uncertainty on A g comparable to that of the 
it^^ir^ decay mode. 

References 



1. J. H. Christenson et al, Phys. Rev. Lett. 13, 138 (1964). 

2. G. Barr et al. (NA31), Phys. Lett. B 317, 233 (1993). 

3. A. Lai et al. (NA48), Eur. Phys. J. C 22, 231 (2001). 

J. R. Batley et al. (NA48), Phys. Lett. B 544, 97 (2002). 

4. A. Alavi-Harati et al. (KTeV), Phys. Rev. D 67, 012005 (2003), Erratum: Phys. Rev. D 
70, 079904 (2004). 

5. B. Aubert et al. (Babar), Phys. Rev. Lett. 87, 091801 (2001). 
K. Abe et al. (Belle), Phys. Rev. Lett. 87, 091802 (2001). 

6. K. Abe et al. (Belle), Phys. Rev. Lett. 93, 021601 (2004). 

B. Aubert et al. (Babar), Phys. Rev. Lett. 93, 131801 (2004). 

7. A. Ceccucci, these proceedings. 

8. L. Maiani and N. Paver, The second DA$NE Physics Handbook, INFN, LNF, Vol 1. 
(1995) 51, 

E.P. Shabalin, |hep-ph/040522"9l Phys. At om. Nucl. 68, 88 (2005), 
A. A. Belkov, A.V. Lanyov and G. B ohm, |hep-ph/0311209| 
I. Scimemi, E. Gamiz and J. Prades, hep-ph/0405204, 
G. Faldt and E. Shabalin, |hep^ph/0503241 , 

9. W. T. Ford et al, Phys. Rev. Lett. 25, 1370 (1970). 

W. S. Choong (HyperCP), Thesis, University of California, Berkeley (2000). 
G. A. Akopdzhanov et al. (TNF-IHEP), |hep^ex /0406008 , 

10. G. D'Ambrosio, G. Isidori and G. Martinelli, Phys. Lett. B 480, 164 (2000), 
E. Shabalin, ITEP-8-98. 

11. R.Batley et al. (NA48/2), Addendum 3 to Proposal P253/CERN/SPSC, CERN-SPSC- 
2000-003, CERN-SPSC-P-253-ADD-3 (2000). 

12. B. Peyaud, 10th Vienna Conference on Instrumentation, Nucl. Instrum. Methods A 535, 
247 (2004). 

13. S. Eidelman et al. (PDG), Phys. Lett. B 592, 1 (2004). 

14. GEANT Description and Simulation Tool, CERN Program Library Long Writeup W5013 
(1994). 



